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We construct the energy-momentum tensor for the gauge fields which describe 



tion laws for energy and momentum, we obtain three different versions for the tensor 
T^'^ , which are physically equivalent. We show that the total energy constructed 



from T^^ is positive for any non-trivial field configuration. Finally, we present a new 
non-abelian solution of the equations of motion for the gauge fields. This solution 
corresponds to spatially uniform color oscillations of the plasma. 
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1 Introduction 

The low-lying excitations of the quark-gluon plasma at high temperature are collective 
excitations which involve the coherent behaviour of a large number of particles over typical 
distances and times of order 1/ gT. These excitations can be conveniently described as 
oscillations of self-consistent average fields to which the plasma particles couple [||]. We 
have shown that, at leading order in the coupling g, the Dyson-Schwinger equations for 
the A^-point functions (A^ > 1) reduce to a set of coupled equations for these mean fields 
and for 2-point functions. The equations for the 2-point functions describe the dynamics 
of the hard plasma particles in the presence of soft background fields; they may be given 
the form of simple kinetic equations for generalized on-shell distribution functions |I]]. 

In the present work, we investigate further the properties of the collective modes, in 
order to properly characterize their energy-momentum content. We restrict ourselves here 
to the case of bosonic collective excitations. Special attention has been payed recently to 
this problem[0, |^, ^. An early attempt to derive the corresponding energy- momentum 
tensor T'^'^ can be found in Ref. 0. There, auxiliary fields are introduced in order to write 
a local effective action, for which the Noether construction is available. The expression 
for T^'^ thus obtained is rather intricate and difficult to use. Another approach is used in 
Ref. , where the plasma is coupled with a weak gravitational field. Again, a complicated 
form is obtained for T^*^, which makes its physical interpretation difficult. This is true, 
in particular, for the terms quadratic in the gauge fields A^, which, however, should 
be analogous to those corresponding to an abelian plasma. We would certainly expect 
the expressions appropriate for the QED plasma to be simple, easy to interpret, and 
closely related to the ones corresponding to classical polarizable media. Another difficulty 
with the derivation proposed in [^ is that the resulting field energy, P*^ = Jd^xT^^, 
appears to be negative even for some simple, plane- wave, field configurations, thus casting 
doubts on the stability of the equilibrium state. It is, however, verified in [Q] that P° is 
positive for the abelian-like plasma normal modes. A similar, but more general, conclusion 
emerges from Ref. [^], where a Hamiltonian for the soft gauge fields is constructed in a 
completely different approach^], by exploiting the analogy between the effective action 
for the collective excitations P, |^, m and the eikonal of a Chern-Simons theory. The 
resulting functional is positive when evaluated for gauge fields which satisfy the non- 
abelian equations of motion without sources, that is, for general, non-abelian, collective 
excitations. 

As we shall show, our approach overcomes most of the difficulties mentioned above 
and has the advantage of offering a transparent physical interpretation. We study a general 



configuration of gauge fields whicli may be induced by an appropriate external color source 
jj^. By using the conservation laws for energy and momentum, we construct the tensor 
j^fiu |-Qj, such gauge fields. In this procedure, we rely on the explicit expressions for the 
retarded non-abelian induced current j*"'^ that we have obtained in previous workslH]. 
The essential aspects of our approach are presented in Section 2. 

Of course, the construction of T^^" from conservation laws does not give a unique 
answer, and, in fact, we shall write down three different expressions, which are physically 
equivalent: 

(i) The first expression, derived in Sec. 3.1. and referred simply as T'^'^ generalizes well- 
known expressions for classical dielectric media (see, e.g., Refs.[p!5|]), to which it reduces in 
the abelian case. The corresponding energy density T°° is manifestly positive for arbitrary 
gauge fields. The same is therefore true for the corresponding field energy P^. 
(ii) The second version for the energy-momentum tensor, referred as T^'^, is derived in 
Sec. 3. 2. It is equivalent, but slightly simpler, to the expression obtained by Weldon in 
Ref . , and differs from it only through a "superpotential" |]l3l . The components T^* 
{i = 1, 2, 3) of this new tensor have a simpler non-local structure than T^^. 
(iii) Finally, the third version, derived in Sec. 3. 3 and denoted as T^^, is traceless and 
symmetric. This allows us, in particular, to construct the gauge field angular momentum. 

In Section 4, we show that the expressions that we have obtained reduce to familiar 
ones when specialized to the case of abelian or quasi-abelian plasmas (i.e., in the weak 
field limit). In Section 5, we briefly discuss new, truly non-abelian, solutions of the 
equations of motion, which correspond to global color oscillations. The main conclusions 
are summarized in Section 6. 



2 Soft fields and hard particles 

We consider an ultrarelativistic quark-gluon plasma close to thermal equilibrum, at a 
temperature T = l/j3 and zero chemical potential. We use natural units, h = c = 1, and 
Minkovski metric. We consider a SU{N) gauge theory with Nf flavors of quarks. The 
generators of the gauge group in different representations are taken to be Hermitian and 
traceless. They are denoted by t" and T°, respectively, for the fundamental and the adjoint 
representation, and are normalized such that tr{t°'t^) = {1/2)6°''' and tr{T°'T'') = N6"'''. 
It follows that (T")bc = —if"''"^, where /"'"^ are the structure constants of the group: 
[T", T''] = if"-^^T'^. We use, without distinction, upper and lower positions for the color 
indices. 



As alluded to in the introductory section, we are studying the response of the plasma 
to a classical, external, color current j^ = j'^T"'- In the absence of the external source, 
the plasma is assumed to be in thermal equilibrium and the expectation values of the 
fields vanish. The interaction with the external current is adiabatically turned on from 
time to ~^ — oo, and the gauge fields acquire then nontrivial expectation values, denoted 
by A^(x). 

We assume the external perturbation to be weak and slowly varying. Then, the 
induced average fields have Fourier components with only soft momenta P ~ gT, and 
their amplitude is constrained so that F^,^ ^ gT"^, or, equivalently, A^ ^ ^Hl- The latter 
limitation on the amplitude of the field oscillations ensures in particular the consistency 
of the soft covariant derivatives: if A^ ~ T , then gA^ ~ gT is of the same order as the 
derivative of a slowly varying quantity, idf^ ~ fi'^^- Occasionally, we shall specialize our 
results to the case of weaker fields, for which it is consistent to replace covariant derivatives 
by ordinary ones; then, A^ behaves as an abelian gauge field and the equations are linear. 
This weak field limit will be referred below as "the abelian regime" . (For a QED plasma, 
the equations remain linear even for fields as strong as allowed, i.e., for A^ ~ T.) 

The gauge fields satisfy the generalized Maxwell equation 

[/^^i^..(a:)r =J»+jmx), (2.1) 

where D^ = d^ + igA^{x), {A^ = A^^T^), and F^, = [D^,D,]/{ig) = F^^T'^. The 
external source j'^{x) is to be considered as a formal device to generate arbitrary gauge 
field configurations. In leading order, the induced current j*"""^ is proportional to the 
fiuctuations of the average phase-space color densities of plasma constituents. These 
fiuctuations are described by color matrices which we denote by 6n±{'k,x) = 6n'^{\<., x) t"" 
and 5A^(k, x) = 6N°'{]i.,x)T°' for quarks, antiquarks and gluons, respectively. Then 



gJ-^^'^{N,[6nl{k,x)-6nl{k,x)] +2N6N%k,x)}, (2.2) 



-—v,TTi^2Nit^[6n+{k,x)-6n4k,x)] + 2T^6N{k,x 



(27r)3 T 

where k^ = \k\ and the factors of 2 account for the spin degrees of freedom. The quantities 
dn± and 6N are determined by the following kinetic equations []T|: 

[v ■ D^, 5n±(k,x)] = T^ V • E(x) ^, (2.3) 

[vD,, 6N{k,x)] = -g^.E{x)^. (2.4) 



Here, v^ = (1, v), where v = h/k is the velocity of the hard particle {k = |k|), E^ = F*" 
is the chromoelectric field. In the r.h.s., Nj. = l/{exp{f3k) — 1) and n^ = l/{exp{j3k) + 1) 
denote, respectively, equilibrium boson and fermion occupation factors. 

Eqs. ( p^) and ( p.4|) generalize the Vlasov equation to non-abelian plasmas. They 
are covariant under a local gauge transformation of the mean fields y4"(x); therefore, the 
fiuctuations (5n^(k, x) and (5A^'^(k, x) transform like vectors in the adjoint representation. 
The gauge symmetry requires the presence of covariant derivatives, which makes these 
equations non-linear with respect to A" If we were to solve these equations for a fixed 
f^ = (1, v), we could get rid of the non-linear terms by choosing the axial gauge v'^A1^{x) = 
0. In this gauge, {v ■ DY^ = S""^ v ■ d and v • E"(x) = — v ■ (SqA" + VAg), as for abelian 
fields. However, in calculating the induced current (|2.2|) we have to integrate over all the 
directions of v. It is therefore necessary to solve eqs. ( p^.3| )-( P~^ in an arbitrary gauge. 
This can be done with the help of the retarded Green's functions for the covariant line 
derivative v ■ D^, defined by 

t{v D^)acG:',t{x,y;v) = 6'^' 6^^\x - y), Gret{x,y;v) = for xq < yo, (2.5) 

and which has the following expression 

Gretix,y;v) = -i9ix'^-y^)5^'^{^-y-vix'^-y'>))uix,y). (2.6) 

Here, U{x, y) is the parallel transporter along the straight line 7 joining x and y, 

U{x,y) =Pexp{-ig f dz'^A^iz)}. (2.7) 

The solution of eq. ( |2.3| ) is then 

d/Tih r^^ 
(5nl(k,x) = Tg—rr duUab(x,x - vu)v ■'E''{x - vu), (2.8) 

dk Jo 

with a similar expression for 6N"'{'k,x). From eqs. ( |2.2| ) and ( |2.(j| ), we readily derive the 
induced current 

/d^ /"oo 
— t;^y^ duUab{x,x-vu)vB\x-vu), (2.9) 

where ujp is the plasma frequency, Up = {2N + Nf)g'^T'^ /18. The integral / dft runs over 
all the directions of the unit vector v. It can be easily verified that jTai^) is covariantly 
conserved, 

D^ f;\x)] = 0. (2.10) 

It follows that the external current is constrained by a similar equation, [D^^, j^lix)] = 0. 



Eq. ( p.9| ) for the induced current summarizes all the information about the collective 
motion of the hard particles in the presence of the soft gauge fields. By using it in the 
r.h.s. of eq. ( ^.1| ), one obtains a generalization of Maxwell's equations in a polarizable 
medium. Besides , the induced current acts as a generating functional for all the (retarded) 
amplitudes between soft gauge fields. For instance, the polarization tensor for soft gluons 
is 

(2.11) 







A=0' 




which gives (with P^ = 


= (/, p)) 




jTab( 


P^-.'^^,2.afe f .0.0 , „0 / 


'dVt 


Vp^Vy 



This is the well-known "hard thermal loop" for the gluon self-energy]^, ^ |1^, [1^. In the 
weak field limit, the tensor II^j, fully characterizes the dielectric properties of the plasma 
For strong fields, A^ ~ T, all A^-gluons amplitudes with N > 2 contribute to the induced 
current (|2.9|) at the same order in gW^ K% M. 



The induced current may be given a different functional form which will be used 
later. We start by defining a new function, W^{x; v), as the solution to 

[vD,,W^{x;v)] = F^P{x)vp, (2.13) 

which vanishes when xq -^ — oo. We have W^ = W^T"", with 



oo 



W^{x;v) = / duUab{x,x-vu)F.^P{x-vu)vp, (2.14) 

Jo 



and the induced current (E.9|) can be rewritten as 



J^aix) =3ujIJ^v,W'^{x-v). (2.15) 

In the case of an abelian plasma, eW^ has a simple physical interpretation: it represents 
the 4-momentum transferred by the electromagnetic field to a particle of charge e and 
velocity v, moving along a straight line from time to —>' — oo to xq. (The transferred 
momentum ~ gT is small compared to the hard momentum /c ~ T, so that the particle 



is not significantly deviated by the Lorentz force.) From eq. (|2.14|) it is obvious that 

v^Wi:ix;v) = 0. (2.16) 

The functions W/^{x] v) are homogeneous of degree zero as functions of v^, i.e., 

.-^'^^^"'=0, (2,17) 



(In taking derivatives with respect to f^, as in eq. ( |2.17| ), we consider the four components 
of f '^ as independent variables, and set v^ = (1, v) only after derivation.) This is obvious 
from the defining equation ( ^.13|) , and can be also verified on eq. (|2.14| ) by using the 
identity 



_d_ 



f{x — vu) = —u{y-dx)f{x — vu). 



(2.18) 



We prove now that the expressions ( p.9| ) or ( p.l5| ) for the induced current are equiv- 
alent to the following one 



3Ta{^) 



82/" dVt 



— ^^ . 
2 n Att 



W::(x;v) +v, 



^v-^i 



dWP{x;v) 
dvP 



(2.19) 



which also appears, with sligthly different notations, in the last of Refs.[0] (as eq. (C.15)). 
To show the equivalence, we note first that eq. ( |2.17| ) implies 



dWP{x; v) 



dW^{x;v) 



'-^^)^°(^'-) 



dvP dv'- 

Then, we use eq. ( p.l6| ) to replace W^{x; f ) by v ■ Wa(x; v) in ( |2.2(]| ), and obtain 



(2.20) 



dW^{x;v) 



W^(x:v) + 



d 



d v'v^ ' 



Wi{x-v) 



The second term in the r.h.s. vanishes upon angular integration, so that 

dn dWP{x; v) r dn 



(2.21) 



An 



dvP 



/£;<(.;..). 



(2.22) 



This proves the equivalence between the time components (/i = 0) of eqs. ( |2.15| ) and 
( |2.19|) . The corresponding proof for the spatial components is similar. 



3 Energy- momentum tensor from conservation laws 

The energy-momentum tensor Tn^{x) of the collective excitation satisfies the following 
conservation law 



d.T^-'ix) = -Fr{x)rJx), 



(3.1) 



for fields satisfying (p.l|) . This equation may be viewed as the thermal average of the 
conservation law formally satisfied by the corresponding field operators in the presence of 
the external perturbation j" (Since no ultraviolet divergence appears at the level of the 



present approximation |p!2|, |I|, we can safely ignore all the complications concerning the 



finiteness of T^^{x) in perturbation theory] 13|.) Note, however, that we shall not derive 



Tnu{x) by directly evaluating the expectation value of some local composite operator 
(e.g., the canonical energy-momentum tensor associated to the QCD lagrangian), but, 
rather, by integrating the conservation law ( p.l|) . This is possible because, in the present 



approximation, all the dynamical information is contained in the equation of motion (^A 
for the gauge fields, and we know, from eqs. ( p. 91) or ( p.l9| ), the explicit expression of the 
induced current in terms of the fields. 

In order to construct the tensor T^'^ , it is convenient to separate the standard Yang- 
Mills contribution, by writing 

T^''{x)^Ti^l,{x) + Q^''{x), (3.2) 

with 

TZA^) ^ -Frix) F:^{x) + \ g^^ Ff{x) F^,{x). (3.3) 

The tensor (|3.3|) satisfies the structural equation 

d.TpiA^) = - [D^ F,,{x) r F7(x). (3.4) 

In deriving this equation, the following identity is useful: 

d^{A^B^) = [D^, A]''B^ + A" [D^, B]\ (3.5) 

We use now the equations of motion ( |2.1| ) for F^'^ to transform eq. ( p.4| ) into 

d.TZ^x) = -F7(x) (j;(x) + j7^"(x)). (3.6) 



By comparing eqs. (|3.1| ), ( p.2|) and (|3.6| ), we conclude that Q^'^{x) satisfies 

a^e'^^(x) = F7(a;)jr'^"(x). (3.7) 



In the next sections we shall use either eq. ( p.l5| ) or eq. ( |2.19| ) for j*"'^ in order to write 



the r.h.s. of eq. (|3.7|) as a total derivative of some gauge-invariant functional of the fields, 
and identify in this way Q^'^{x). 

The tensor T^^, thus defined will satisfy, by construction, the conservation law ( |3.1| ), 
and will vanish for large negative times. Accordingly, the 4-momentum constructed from 
it. 



P''{t)= fd^xT^^it,^), {3.i 



will coincide, at any time, with the correct 4-nionientuni of the collective excitation, which 
is obtained by integrating the conservation law 

doP^t) = ~ Jd'xFr{x)j;{x). (3.9) 

(In priciple, this equation determines P'^{t), given the initial condition P^(t —>■ — cxd) = 0. 
Recall that the external sources and the average fields are assumed to vanish at t ^ — oo). 
Of course, the tensor T^,^ thus obtained is not unique, since any two tensors which would 
differ by a divergence-free piece, and lead to the same P^, would be equally acceptable. 
We shall indeed present in this section three different versions of the energy-momentum 
tensor. 



3.1 A positive definite energy density 

Eq. (|2.15|) for j™'^ allows us to write the r.h.s. of eq. ( p.7|) as 

= ^^ij^ [vD^, W^{x-v)r W',{x-v), (3.10) 

where eq. ( |2.13| ) for W^x^v) has been used in the second line. For z/ = 0, this can be 
easily written as a total derivative. 



W',{x;v) 



^^{x)--C\x) = 3ulJ^[vD^,W\x;v) 

= d.e^'^ix), (3.11) 



with the definition 

dn 



Q''{^) ^ ^u;lJ — v^W'^{x;v)W'^{x;v), (3.12) 

where the identity ( p.5| ) was used. After adding Tym{x), eq. (p.3|), we end up with the 
following expressions for the energy density 



T^"" x) 



^(E'^(a:)-E^(x) + B^{x)-B%x)) + ^^j/ ^ C(^; t.) C(x; .;), (3.13) 



and for the energy fiux density, or Poynting vector, S^ = T 

39/" dQ 



iO 



S{x) = E%x) X B^{x) + ^^lj^^ W'^{x;v)W'^{x;v). (3.14) 

The chromomagnetic field components are, as usual, B\{x) = —{l/2)e^^^Fi^{x). 



The expression ( p.l3| ) for T°°(x) is manifestly positive semidefinite and the same is 
obviously true, at any time t, for the excitation energy S(t), 



Sit) = j (fxT'^'^it, X) = SYM{t) + Wit). 



(3.15) 



Note that 

£YM{t) = Jd^'^l (E"(a;) ■ E%x) + B%x) ■ B%x)) (3.16) 

also involves medium effects, since E"(x) and B"(a;) are the total gauge fields, including 



the plasma polarization. The second piece of (3.15) 



W{t) ^\ujI j d'x j '^^W!{x,v)W'^{x,v), 



(3.17) 



may be interpreted as the polarization energy oi the plasma, that is, the energy transferred 
by the gauge fields to the plasma constituents as mechanical work of the chromoelectric 



field. This becomes more obvious if eq. ( p.ll|) is used to reexpress W(t) as 



W(t) 



dt' / d^xW{t\ x)-j™'^(t', x). 



(3.18) 



One can easily verify - by explicitly performing the integral over t' in the r.h.s. - that 
eq. ( p.l8|) leads indeed to the expression (|3.17| ) for Wit) if eq. (|2.15| ) for the induced 
current is used. 



The positivity of the polarization energy ( |3.17| ), and hence of the total excitation 
energy (|3.15|) , reflects the stability of the equilibrium state with respect to long wavelength 
color fluctuations. 

We briefly discuss now the other components of the energy-momentum tensor G^'^(a;). 
We can write eq. (p.7|) as 



d.Q'^'ix) 



dVL 



F:^ix)J",ix-v) 



(3.19) 



where, according to eq. (|2.15|) , J^ix;v) = ?>uj'^v^W^ix]v). We denote by Wix^v) the 
solution to 



ivd.,)Wix-v) = F7(a;)J-;(a;;t;), 
which vanishes for xq — > — oo, that is. 



(3.20) 



Wix-v) 



t/rF7(x - vt) J^ix - VT- v). 



(3.21) 



It is then easy to see that 

e'^^(x) = f^v''U''{x;v) (3.22) 

satisfies eq. (|3.19| ). For z/ = 0, the integral over r in ( p.21|) can be easily performed, with 



the result n°(x; v) = (3/2) u^ Vr°(x; v) Wj^^x; v). By inserting this result in eq. (|3.22| ), we 
recover our previous expression for 0'^°(x), eq. ( p.l2| ). The physical content of the function 
Il'^{x] v) is transparent: it represents the density of the total 4-monientuni transferred from 
the fields to the hard particles with velocity v from time — ;► — oo to Xq. In particular, 
B°^(x) is simply the angular average of 11^ (x; v) over the possible orientations of v. 

The expression ( 3.22| ) for B'^''(x), though formally simple, involves generally a double 



integral along the hard particle trajectory (see eqs. ( ^.21 ) and (2.14)). We shall construct 



in the following section a different expression for the energy-momentum tensor, which 
involves a single integral along the trajectory. 

3.2 Another expression for T^^ 

We use now the second version for the induced color current j"^^, eq. ( p. 191 ), in order to 
construct another energy- momentum tensor, T^^'^{x) = Ty^Mi^) + ©'^'^(a;), different from 
the one obtained in Sec. 3.1, but physically equivalent. 

When using eq. (|2.19| ) for jj^^, the r.h.s. of eq. ( |3.7|) reads 
F7{x)f;'''{x) = 1^:1 J ^ 1^7(0:) W^^{x; v) + F7(a;) v, ^^'J^J; ^^ | . (3.23) 

In the first term inside the brakets, we express F^^{x) in terms of the functions M/^(x; f ), 
using the identity 

[D^ W^'i.x-.v)] - [L"^, Wi^x-v)] = Ff'^ix) + ig[W^'{x;v),W''{x;v)]. (3.24) 

(To prove this identity, act on eq. ( |2.13| ) with the covariant derivative D'^. Then, consider 
the similar equation where the indices fi and u have been interchanged, and take the 
difference of the two equations. The identity (|3.24| ) follows after some repeated use of the 
Jacobi identity.) We obtain 

Fn^) W;ix; v) = [D'^, W%x- v)^ W;{x- v) - ]^d^ {W ■ W) , (3.25) 

with the notation W -W = Wj^lx; v)W^{x; v). For the second term in the r.h.s. of (|3.23|) , 
we use eq. ( p. 131) for W"^{x] v) to write 



dWP dWP 

10 



/ dWP\ 



dvP ) 



v-D, 



dWP 
dvP 



By using eq. (|2.13|) once again, we see that 



d 

dvP 



[vD,,WP{x;v)Y 



d 

dvP 



{Fn^)v, 



so that 



dWP 



dWP\ 



" ^ dvP ^'^ V" "" 9vP J 

By adding together eqs. ( p.25| ) and ( |3.28| ), we finally obtain 



Fri^)j, 



indai 



[X 



d^Q^''{x) 



with 



e^^fx) 



dVL 



dwi 

" dvP 



1 
2 



The polarization energy density associated to ( |3.3U| ), which is 



(3.26) 



(3.27) 



(3.28) 



(3.29) 



(3.30) 



(3.31) 



is different from the previous expression in Sec. 4, eq. ( p.l2| ). In particular, the expression 
above is not manifestly positive. One advantage of the new tensor above is that eq. (|3.30| ) 
offers, for the components (//,«), a simpler non-local structure then eq. ( |3.22| ). 

It is possible to show that the energy- momentum tensor TP'^{x) is physically equiv- 
alent to that proposed by Weldon in Ref.0. Indeed, the difference between the two 
tensors is the total derivative of an antisymmetric tensor of rank 3. To be specific, we 



haver'^''(x)-T^^(x) 



dpXPP^{x), with 



XPP-'ix) = -uo^ 



dVL 



[(vPV^,{x-v) - v^VP,{x-v))w:{x-v)}, 



(3.32) 



2 " P 7 47r 

and XPP^{x) = —Xpp^{x). Here, TiJ^{x) denotes Weldon's tensor and the functions 
V^(x; v) satisfy, by definition, the following equation: 

[vD^.V^ix-v)]" = W^{x-v). (3.33) 

Eq. ( p.32|) is obtained by first transforming Weldon's original expression with the help of 
the identity (|]2|). 

The tensors Qp^ and O'^'^ constructed up to now are not symmetric, which is not 



surprising since the indices jjl and v enter disymmetrically in the conservation law (|3.7|) 
on which our derivation is based. 
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3.3 A traceless and symmetric tensor 



We construct now a symmetric tensor, starting again from eq. ( |3.7| ), but using in its r.h.s. 
the expression (|2.2|) for the induced current, which we write as 



jr''°(^) 



«/7m^^"(^-) 



with 



Sr{k,x) = Ni[Snl{k,x) - 6n''_{k,x)\ + 2NSN''{k,x 
The r.h.s. of eq. (|3^ ) reads then 

Frix)f;"'^{x) = g I ^—^^F:'^{x)v^6nk,x). 



i2ny 



(3.34) 



(3.35) 



(3.36) 



We use the simple identity 



F:^ix)v, = F:^{x)^^{v'^k,), 



(3.37) 



(recall that v^ = k" /k^ and k^ = |k|), and an integration by parts, to rewrite ( p.36|) as 



Fi:nx)f:^^{x) = -g I ^^^^v'FT{x)k^^6r{\i,x). (3.38) 



a \ / J fi 



{2ny 



d 

"^ dk^ 



Let Se(k, x) be the solution to 



{vd,)6e{k,x) = -gFT{x)k, — 5nk,x) 



d_ 

dk^ 



which vanishes when xq ^ — oo. Then 



5e(k,x) = iVf 



(5e+(k,a;) + (5e_(k, x) + 2N 5e^{k,x) 



with obvious notations, and eq. ( ^.38] ) becomes 

F7{x)3^;^'^{x) = 9,e^;^(x) 

where 

d^k 



(3.39) 



(3.40) 



(3.41) 



e^;^(a;) = / j^^v^v-'{Ni[5e+ik,x) + 6e4k,x)] + 2N5e,ik,x)} (3.42) 



i2ny 



is a traceless symmetric tensor. The same is therefore true for the corresponding energy- 
momentum tensor T^^^^x) = Ty^j(x) + Q^^^^x). Remark that, in the weak field limit, 
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Se(k,x) is quadratic with respect to the gauge fields (see eq. ( |3.39| )), as opposed to 
(5/"(k, x) which, in the same hmit, is hnear in A^. 

The symmetry oiT^^{x) allows us to construct a simple density for the generalized 
angular momentum of the gauge fields: 

M"^%x) ^ x^ T^;^{x) - x^ T^/^ix). (3.43) 

It can be verified, by using eq. ( |3.1| ), that M'^'"^(x) satisfies the correct conservation law 

dxM^'^'ix) = - (x^F;^(x) - x''F^\x)) j^ix), (3.44) 

where one recognizes in the r.h.s. the torque applied to the system by the external current. 
Since M^^^(a;) vanishes when xq —>■ — oo, the total angular momentum is correctly given 
by 

JA'^(t) = fd^xM^^'^ix), (3.45) 

for any t. 

The vanishing trace of T^^^^lx) refiects the dilatation invariance of the massless 
tree-level QCD lagrangian. At the level of the present approximation, there is no break- 
ing of the dilatation symmetry, neither by radiative corrections (via the regularization 
procedure), nor by a thermally induced mass, which is negligible for the hard particles. 

4 Applications 

In this section we use the general expressions derived in Sec. 3 to study the energy transfer 
in the plasma and the energy-momentum content of particular excited states. We consider 
mostly the weak field, or abelian, limit, where the components of T^'" are quadratic in the 
gauge potentials. We show that, in this particular limit, our expressions are the straight- 
forward generalization of the corresponding formulae in the standard plasma physics. 

Moreover, any solution of the abelian gauge field and kinetic equations can be 
directly imbedded in the corresponding non-abelian equations. That is, let A^{x) and 
W'^{x;v) be solutions to the coupled abelian equations 

d^F,,ix)-3u;lJ^v,W%x;v) = j,(x), (4.1) 

and 

{vd^)W\x;v) = v-E(x). (4.2) 
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Then, A^(x) = X°- A^{x) and W^{x;v) = X"-W'^{x;v), with constant A", satisfy the cor- 
responding non-abehan equations, (|2.1|) and ( p.l3|) , for an external source pointing in a 
fixed direction in color space, jf^ix) = X°'j^{x). This property extends the corresponding 
one for classical Yang-Mills equations in vacuum (see, e.g., Ref. |jl4 



A slightly more general field configuration for which the general solution is still 
trivial corresponds to gauge fields in the Cartan subalgebra of SU{N), that is, A*(t) = 
J2s ^l{t)T^, where the indice s takes on A^— 1 values corresponding to the A^— 1 commuting 
generators T" of SU{N). Then, clearly, F^^ = d^^A^ - 9^v4^, and F^" = W^ = j^ = 
for a y^ s. The non-trivial components A^, W^ and j^ satisfy abelian-like equations, as 
eqs. ( |4.1| )-( |^?^ ) above. Hence, the problem reduces itself to A^ — 1 independent abelian 
theories. 

In subsection 4.1, we briefly discuss the case of static gauge fields. In the following 
subsection, we investigate the energy transfer between fields and particles and, in partic- 
ular, the energy dissipation through the Landau mechanism. Finally, in subsection 4.3, 
we consider some interesting special field configurations, namely, time-periodic fields and 
abelian plasma waves. 

4.1 Static field configurations 

We consider a field configuration induced by a static external color density p"(x). We 
choose a gauge where the potentials A^ are time-independent. In this case, it is useful to 
reexpress W^{x; v), eq. ( |2.14| ), as a function of the gauge potentials themselves 

W°{x;v) = -A'^{x) + duUab{x,x-vu)v- A\x-vu), (4.3) 

Jo 

where A" (a;) = (9ov4"(x). For a static field configuration, A"'{x) = and 

C(x;^) = -A°(x). (4.4) 

It follows that the induced current ( p.l5| ) reduces to a static color distribution 

j;"'(x) = - ^,0 3 ul A°(x) = - g,o P':\x). (4.5) 



Then, the gauge field equations ( p.l|) become simply 



[A,i?^(x)] +3^,^A°(x) =p(x), 

e''\Dj,B^] = ig[A'',E']. (4.6) 
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The first equation above clearly shows the Debye screening of the static chromoelectric 
field, with the screening length Xd = 1/(3 cUp)^/^. Note also that the chromomagnetic field 
induced by a static color charge is generally not zero in a non-abelian theory. Furthermore, 



n -^p - -a\--/ - -a\--/ r\ 



e°°(x) = -^P^a(x)A°(x) = --pr'(x)A°(x). (4.7) 



In analogy with what is commonly done in abelian plasmas, (see Appendix), we define 
the color polarization vector P"(x) by [D*,P*(x)]'^ = p^"'^(x), and construct the color 
displacement vector 0*^(0;) = E"(a;) + P'^(x). Then, the polarization energy (|3.17| ) can be 
written as 

W = fd^xe^\^) = -- I d^x[D\P\-K)YAl{-K) = - /'rf=^xE'^(x) ■P'^(x), (4.8) 

as familiar in electrostatics of classical dielectric media[]T5|. (We used here -E*(x) = 
[D^, A^lx)]"- for static fields, together with eq. (|3.5| ) and an integration by parts.) The 
total field energy follows: 

S = - fd^x (E"(x) • D"(x) + B'^(x) ■ B"(x)). (4.9) 

For an abelian plasma, we can solve the equation 



to get the gauge field 



-A + 3cu2)A0(x) = p(x) (4.10) 



A°(x) = f^PMe-mn^ (4.n) 



where R = |x — y|. Then, the field energy can be expressed in terms of the external 
charge, with the standard result 



2 J (27r)3 p2 + 3^2 • 



S = -Jd'x E^(x) + ^ulAl{^)) = ; /t^ T, o.. - (4-12) 



For a non-abelian plasma, the non-linear structure of the generalized Maxwell equa- 
tions (|4.6| ) does not allow such a simple elimination of the gauge fields in terms of p". 



4.2 Energy transfer and Landau damping 

When studying time-dependent field configurations, there is a new feature arising: the 
imaginary part of the polarization tensor ( |2.12| ) is non-vanishing, which indicates that 
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there is dissipation of the gauge field energy in the medium. To investigate this, both for 
abehan and non-abehan fields, we compute the rate of energy absorbtion by the plasma 
constituents (see eq. ( |3.18| )): 

dW{t) 



dt 



Jd'xE^{t,^)-ir%^). 



(4.13) 



Consider first an abelian-like plasma, in which case the general description in Sec. 2 
reduces to the standard Maxwell equation in a polarizable medium||16[|. For simplicity, 
consider a time-periodic field 



E(t, p) = E(p)e''* cosu;ot 



(4.14) 



with 7/^0+ and both ujq and p = |p| of order gT. Then, a simple calculation shows that 
the average energy loss of the field, that is, the expectation value of eq. ([4.13|) over the 
period Tq = 27t/ujo, is proportional to the imaginary part of the dielectric ten5or|T6|: 



'dW' 
~dt ^ 



d^p 

f2^ 



^oi^'(-p)(lme^^(a;o, P) E^ {p) 



(4.15) 



(The definition of e^\ as well as its relation with the polarization tensor, are reviewed in 
Appendix). In our case, eqs. (|A.12|) and ( p.l2| ) show that 



Ime*-'(c(j, p 

and we obtain the average energy loss as 
/dW\ 



Sir —^ — V v' oiio — v ■ ij), 
uj J Att 



dt 



3 2 



/7^/S^(--p)K-(p) 



(4.16) 



(4.17) 



This expression is obviously positive: on the average, the fields transfer energy to the 
particles. This energy dissipation is related to the Landau damping of the mean fields]!^, 
that is, to the coherent transfer of energy in between the gauge fields and the hard particles 
which are moving in phase with the field oscillations. As the 6 function clearly shows, the 
dissipation is kinematically allowed only for fields carying space-like momenta (i.e., when 
Uo = v-p). 

The energy given by the fields to the colored particles is ultimately supplied by the 
external source, that is 



'dW' 
dt , 



/ d^xE{x) ■ i{x] 



(4.18) 
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as shown by the conservation law ( |3.1| ) (the local piece dSyM/dt averages to 0). Therefore, 
in order to compute the energy loss of a given external source, one just has to replace the 
electric field E"^ in eq. ( 4.15 ) in terms of the current j^, by solving the Maxwell equation 



( |2.1| ). In the abelian regime, this can be easily done by Fourier analysis, and leads to 
the following expression for the average energy loss of the time-periodic external current 
j(p) cosujot: 

'dW\ 1 / d^p [ |p(p)|\ 1 , |j.(p)|\ 1 1 ,..., 

-^oi 3 — ^^—^ tH 2 — Im— -—^^.(4.19) 



\ dt / 2 J (27r)3 I p2 ei{uJo,p) Uq et{uJo,p) - p'^/ujq 

Here, uq p(p) = — «p-j(p) , as required by the current conservation, while jt(p) is the trans- 
verse component of j(p), i.e., jj(p) — i^^'' ~ V^V^) j^iv)- The longitudinal and transverse 
pieces of e'-' are defined in Appendix (eqs. ( |A.13| )- ([A.15| )). 

Clearly, eq. (|4.19|) also holds for non-abelian field configurations such that the space- 



time and the color degrees of freedom decouple, that is, for the imbedded abelian solutions 
to eq. ( p.l|) with a color source j'^(t,p) = A'^j(t, p), and constant A". (The only modifi- 
cation consists in a supplimentary factor A"A° in the previous equation.) An expression 
similar to ( [4. 191 ) has been used in Refs. ]TB| , |I^ in order to compute the energy loss in 



the plasma of an external parton, which is assimilated to a classical point-like colored 
particle, moving with constant velocity v, and therefore giving rise to a current 

■f{uj, p) = 27r A" V (5(tu - V • p). (4.20) 

Clearly, the energy loss calculated according to eq. ( |4.19|) accounts only for the energy 
transferred by the external source to the hard particles, through its coupling to the collec- 
tive modes. This is indeed the dominant process for long wavelength and low frequency 
excitations. However, these conditions are not satisfied for the whole range of momenta 
involved in the source (|4.2CI| ). Strictly speaking, when used with such external sources, all 



our previous integrals over p should involve an ultraviolet cut-off of order gT. 

To close this subsection, we return to the general case of an arbitrary non-abelian 
gauge field configuration (with \A°j\ ^ T) and compute the total energy absorbed by the 
plasma during the life-time of the excited state: 



AW 



dt-P- = dt d^xB^it, x) ■jt'^it, x). (4.21) 

-oo d t J —oo J 



This quantity is well defined if we assume that the external perturbation is switched off 
adiabatically at t -^ oo. It coincides, of course, with the t — > oo limit of the polarization 
energy W(t) given by eq. ( |3.18D . By using eq. ( |2.9| ) for j^"'^, we obtain, after some 
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elementary changes in the integration variables 
3 o /■ dQ 



AW = -cjM ^ / rf^xv-E"(x) / duUab{x,x-vu)w^\x-vu). (4.22) 
2 J An J J-oo 

This quantity vanishes if 

duUab{x,x — vu) El{x — vu) = 0, (4.23) 



for arbitrary x, v and i. For fixed v, we may choose the axial gauge v'^A"^ = 0. Then 
eq. (|4.23|) becomes 







duElix — vu) 



(2^ 



-^P-^ArnO 



5(p'^-vp)E:(/,p). 



(4.24) 



As it should be valid for any x, eq. ( |4.24|) implies El{p^ = v ■ p, p) =0, which is the 
condition for absence of Landau damping. More generally, eq. ( [4.23|) characterizes the 
non-abelian gauge fields which propagate without dissipation. It is one of the conditions 
invoked in Ref . M in order to write a well-defined effective action for the soft gauge fields. 



4.3 Plasma waves 



We start by reconsidering the periodic field ( [4.14|) and assume that E(p) = if |p| > c^o, 
so that there is no dissipation. The average, over the period To, of the polarization energy 
Wit), eq. (^), is 

2 



(W) 



TCUJ^ 



d^p fdn V ■ E(p) 



(27r)3 J Ait {loq - v ■ p) 



The angular integral can be easily evaluated (see eqs. ( |A.21| ) and ( [A.l 
the Yang-Mills contribution. 



i^YAl) 



d^p 

(2^ 



E(P) 



+ 



UJn 



P X E(p) 



(4.25) 
After adding 

(4.26) 



we derive the following expression for the averaged value (over Tq) of the total field energy 

s{t) = SYM{t) + w{ty. 



{£) 



d^p 
A J {2nf 



OJ, 



l + 3^a(a;o/p) 



ekp: 



+ 



2 2 

1 + ^ + 3 ^ K^o/P) 
^0 V 



E,(p) .(4.27) 



The coefficients aiio^jp) and hiuj^/p) are given in eq. (|A.23|) . They are both positive for 
Uq/p > 1. Furthermore, E(p) = Ej(p) + Ei(p), with p ■ Ei(p) = 0. 
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The coefficients in eq. ( |4.27] ) can be related to the derivatives of the dielectric ten- 
sor (|A.12|) with respect to u. Indeed, by using eqs. ( |A.14| )- ([A.15| ) for ei^tiy^^p), it is 
straightforward to verify that eq. ( [4.27[ ) coincides with 



<^)4 



d 



(27r)3 \ dujQ 



uJoei{uJo,p) Ei{p 



+ 



d 

dujQ 



uJo{et{uJo,p)-p'/uj'Q) Ei(p) .(4.28 



A similar expression holds for ordinary polarizable niedia[|I^. 

The previous discussion can be easily generalized to study the abelian plasma waves, 



i.e., the solutions to eq. (|4.1| ) without external sources. The decomposition of an arbitrary 
electric field in terms of normal modes is of the form 



E(t, x) = /yl^e^P- Y: {e.(p)e"^-(^)* + e 



s=Lt 



-p) e'^'^P^'} 



(4.29) 



where the frequencies ujs{p) are solutions of the dispersion equations 
ei{uJi{p),p) =0, and et{ujt{p),p) = p'^/uj'^, 



(4.30) 



for longitudinal and transverse modes, respectively!^, [T^, |2^. Note that 07/(0) = u!t{0) = 
Up and uJi^tip) > P, for any p, and hence the normal modes propagate without damping. 

By using the dispersion relations ( |4.30| ), and eqs. ( |A.17| )-( [A.18| ), we easily obtain 
the energy of the field ([4.29|) as 

d^p [ / 



S 



H 



(27r)3[V^F(p)-p^ 



+ 



^ujI 



+ 



p 



lleKp)-e;(p) 

,2 \ 



1 e,(p) 



IP) 



(4.31) 



,^t{p)-p'^ ^Kp) 

As expected, this is time-independent. The coefficients in eq. ( [4.31| ) are inversely propor- 



tional to the residues z^tip) of the effective propagators ( |A.4|) at the appropriate poles 
(see eqs. ( [A.6[ ) and ( |A.20D ). Since the energy ( |4.31| ) is positive, the same holds for the 
residues zi^tij>)] this can be also verified directly by using the formulae given in Appendix. 

It is quite straightforward to verify that the second expression for the energy density, 
i.e. eq. ( ^.31| ), leads to the same excitation energy. More interestingly, eq. ( |3.30| ) allows 
us to compute easily the total 3- momentum of the field configuration ( [4.29|) . We obtain: 



P 



?>^l 



d^p 
(27r)3\^/(p) \ujf{p)-r 



-1 eKp)-e;(p) 



P 



3cj2 



+ 



p^ 



^t{p) \ujt{p) -p^ ^Kp) 



1 e,(p) 



IP) 



(4.32) 
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We have used here eq. ( [A.22| ) in Appendix, as well as the dispersion relations for the 
normal modes, eqs. (^4.30|) . 

Physically more transparent expressions for the energy and the momentum of the 
plasma waves are obtained by expressing them in terms of the gauge potentials which 
correspond to the field ([4.29|) . We use the normal modes decomposition 



^'^^^^ = / tC^ ^ 4=¥v {^'^(P' ^^ «a(p) e-^^- + e; (p; A) a*, (p) e^^^} , (4.33) 
•^ \^^) A=o y2ujx{p) 

where A = corresponds to the longitudinal mode, A = 1, 2 to the transverse ones, 
p^ = ujx{p), and the polarization vectors e^(p; A) are defined in Appendix. In terms of 
ax{p) and ^^(p)' ^^^ plasma wave energy and momentum become simply 

^ = /(0i|^Kp)«Kp)ar(p)+^t(p) E^«A(p)al(p)|, (4.34) 

and, respectively, 

P = /(0^{p«KP)«r(p) + PE^aA(p)aI(p)|. (4.35) 

The interpretation of the two equations above in terms of the color elementary oscillations 
is obvious. 



5 Global color oscillations 

In this section we study particular non-abelian solutions to the field equations which 
are uniform in space, i.e., which describe global color oscillations of the plasma. For 
convenience, we choose the temporal gauge, v4°(x) = 0. Therefore, v4^(a;) = (0, A"(t)). 
For such fields, it is useful to work with the following representation for lV°(x; v): 

W^{x;v) = -Y- A''{x) - r duUab{x,x-vu){Y-V^)(y- A\x-vu)). (5.1) 

This can be derived from eq. (|4.3|) , by noting that U{x,x — vu) satisfies the equation 

d 

—- U(x,x — vu) = —igU(x,x — vu)v-A(x — vu), (5.2) 

ou 

where A^^ = A'^T"'. For fields which do not depend on spatial coordinates, we have simply 

W^{t;v) = -vA'^(t). (5.3) 
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From eq. ( |2.15| ), one obtains then p^^'^{t) = and Ja"'^(t) = — cj^ A"(t), and, using the 
definition ( gA2D , 



e'^it) 



-ulA.\t)-A\t). 



(5.4) 



This expression of 0*^° is very similar to that obtained in the static case, eq. ( [4.7|) ; in 
particular, both eqs. ( |4.7| ) and ( |5.4| ) are local functionals of the fields. The color electric 
and magnetic fields are 



E'^(t) 






B^it) = -r^''A\t) xA%t). 



(5.5) 



Using these expressions, we obtain a simple expression for the total energy density in 
terms of the vector potentials 

TOO(t) = If^.^ + ^^^a'^-aA + ^r'T''\A''-A''){A''-A'). (5.6) 

Eq. (|5.4| ), together with the other components of the energy-momentum tensor, can also 
be obtained from O^y^ given in Sec. 3. 3. Indeed, for A'^^x) = (0, A"(t)), eqs. (p.39|) and 
( Wm imply 



/^,A„^{2(v.A")(v.A")-iA".A»} 



sym PI 4^ 

and it is easily verified that 0°™ = 0°°. Furthermore, 0g*„ = and 



(5.7) 



sym 



0'J 



sym 



I 4 ^^^ I 

coll-Al^Ai-^A'^.A^. 
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(5.^ 



The external color sources which are needed to generate the field configuration can 
be determined from the equations (|2.1|) which take here the form 



and 



fit) 



Pit) 



cPA' 
IF 



^9 



A' 



,. dA 



dt 



+ ^;A' + g' 



A\ A 



A 



(5.9) 



(5.10) 



Note that these equations differ from the corresponding classical Yang-Mills equations in 
the vacuum only by the presence of a mass term ~ iJ^ A^ for the vector potential. 

In what follows, we look for self-sustained global color oscillations, i.e., solutions 
to the equations (|5.9|)-([^?T(]|) with vanishing external sources (p(t) = fit) = 0). The 
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corresponding equations in the vacuum have been extensively investigated, mostly for the 



SU{2) color group|21|. However, as we shall see, the presence of the mass term oOp ~ gT 
dramatically affects the dynamical content of these equations. Besides, the respective 
solutions for the high temperature quark-gluon plasma have direct physical relevance, as 
they represent possible collective excitations of the system. The physical interpretation 
of the corresponding excitations in the confining regime is more obscure. 

The condition p(t) = implies, according to eq. ( |5.9| ), the commutativity of the 
color matrices A^{t) = A^(t)T" and v4*(t) = {dA\/ dt)T"' . We shall consider here two 
non-trivial possibilities: 

(i) The vectors {^^(t)} and {A\{t)} are parallel in color space for any i, which implies 
A\{t) = Al^h^{t) (no summation over i), with constant Al^ and arbitrary functions /i*(t). 
(ii) The color matrices A' (t) belong to the Cartan subalgebra of SU{N) , A* (t) = J2s Al {t)T^ 
(Recall the discussion at the beginning of Sec.4.1.) If A^ = 2, (ii) is just a particular case 
of(i). 

We consider now the consequences of the second set of equations, that is, j*(t) = 
for all i, on the two possibilities mentioned above. 

(i) For gauge potentials of the form A^{t) = A^K'it), [A^ = A^T"'), the normal 
modes equations f = give generally non-linear and coupled second-order differential 
equations. As an example, consider SU{2), where /"'"^ = e"^'^, (a = 1, 2, 3), and the 
particular field configuration in which A^^ = 5^. Then, the functions /ij(t) satisfy 

^+ujlh+g'hihl + hl) = 0, (5.11) 

plus two similar equations for h2(t) and hs(t). The associated energy density, 

T°° = ^e[(^)' + ^>.^] +y(/^?/^2 + /^?/^^+/^^/^9, (5.12) 

is an integral of the motion and acts as an effective Hamiltonian for the functions hi{t). 
Contrary to what happens at zero temperature, where Up = 0, here the energy conser- 
vation prevents any trajectory {hi(t)} from getting too far away from the origin in the 
3-dimensional space with axis hi. 

Periodic solutions to eqs. ( |5.11|) can be easily written down in the particular case 
where all the three colors oscillate in phase: hi = h2 = h^ = h ("white oscillations"). The 
function h{t) satisfies then the non-linear equation 

h + culh +2g^h^ = 0, (5.13) 
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which has the following solution 

h{t) = hecn{neit-to);k), (5.14) 

for the initial conditions h{to) = hg and /i(to) = (the overdots indicate time derivatives). 
Here, cn(x; k) is the Jacobi elliptic cosine of argument x and modulus /c, and to is the 
arbitrary origin of the time. The quantities k, hg and VLq are related to the dimensionless 
parameter 9"^ = [g"^ /uoT)T^^ by 



k = ^ 



3 



a;„ 



and 



V2^ 



fifi 



CJ, 



(1 



The solution (|5.14|) is periodic, with period 

4 



% 



fifl 



K{k), 



1/2 



,2^V4 



nl/2 



(5.16) 



(5.17) 



(5.1^ 



where Kili) is the complete elliptic integral of modulus k. Since \h{f)\ ^ T, then 6* ^ 1 
and Tg is of order of 7J) = 271 /ojp. The associated field strenghts are ii^* = 
Bl^ = gh? 61^. Therefore, the vectors E^ and B^ are parallel for any a. 



61^ h, and 



Other non-abelian solutions to eqs. (|5.11|) are discussed in Ref. ||2z 

(ii) For gauge fields in the Cartan algebra, A\t) = J2sAl{t)T'^: the commutator 
entering eq. ( |5.1U| ) is trivially zero, so that the normal mode equation reduces to 



df" 



+ ^;K 



Its general solution has the form 

Kit) - 



Al cos UJpt 



0. 



Bl sin UJpt, 



(5.19) 



(5.20) 



with arbitrary constants Al and Bl- A particular solution of this type for SU{3) is 
A(t) = T^ acosujpt + T^hsinujpt, where a and b are fixed vectors in coordinate space. 
This describes coupled oscillations in both coordinate and color spaces, with an energy 
density 

1 



T°° = ^ul{si' + h'). (5.21) 

In particular, a* = 6* corresponds to global color rotations in the plane T^ ~ T^. 
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6 Conclusions 

In this work we have studied the energy and the momentum of a general soft gauge field 
configuration describing collective excitations of the quark-gluon plasma. Our derivation 
relies on the explicit formulae that we have obtained previously for the induced current 
expressing the plasma response to the soft gauge fields. The simplicity of these expressions 
allowed us to construct three different, but physically equivalent, versions for the tensor 
T^^", which are gauge-invariant and rather simple. One of these versions is symmetric (and 
traceless), which allowed us to simply obtain the gauge field angular momentum. We have 
shown that the total excitation energy is positive, for arbitrary gauge fields. We have also 
verified that, for particular abelian field configurations, our general expressions reduce to 
well-known formulae from classical plasma physics. Finally, we have analyzed the global 
color dynamics of the plasma and given a particular non-linear solution which describes a 
truly non-abelian excitation (the plasma white oscillations). The specific examples that 
we have considered in the last two sections illustrate, in particular, the utility of our 
general expressions for T^^" for various applications. 



Appendix 

In this appendix we discuss the leading order effective propagator for soft gluons, 
together with some related topics, like the dielectric permitivity tensor for gauge plasmas 
and the polarization vectors for the gluonic normal modes. 

As already discussed in Sec. 2, the dominant contribution to the polarization ten- 
sor for soft gauge fields is given by Il^i^{P), eq. ( p.l2| ), and satisfies P'^U^u^P) = 0, 
consequence of the conservation law (|2.10|) . The soft gluon propagator is therefore 
*D~l{P) = DqI^{P) + n^i.(P) fTT], [T^. In order to invert this expression, it is use- 
ful to decompose 11^,^ (P) into longitudinal and transverse components (relative to p)|TO[]: 

li^'^iP) ^ Mp^pwr + nt(/,p)pr, (A.i) 

where the projection operators Pf^ satisfy 

p^vr = p^.vi^'' = 0, vr + vr = (p^pyp^) - g^", 

vr Vip. = - vt,, vr Vtpu = - vr, vr vtp. = o, (a.2) 

with P"^ = Pq — p^ and p* = p'^/p. The explicit forms of the functions Ili^t{p^,p) follow 
easily from eqs. ( [A.1| )-( [AI2D and ( |2.12| ); they can be inferred from eqs. ( [A.14| )-( |A.15| ) 
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below. Then, the gluon propagator is 

*D/^-(p) = -vr *MP) - vr *MP) + \ ^^, (A.3) 

in a covariant gauge with gauge fixing parameter A (A = 1 in Feynman gauge). Here 

*Ar\P) ^p' - nK/,p), *A-i(P) ^ p2 - n,(p°,p), (A.4) 

are the effective inverse propagators for longitudinal and, respectively, transverse gluons. 
They vanish on the mass-shell for soft quasi-gluons, thus defining the corresponding dis- 
persion relations, which we denote by ±u;s(p), with s = I or t. The mass-shell residues 
are defined by the relation 

*^siuj,p) ^ ^^— ^ ^"'^^l (A.5) 

where the approximate equality holds in the vicinity of the pole. It follows that 

As discussed in Sec. 4. 3, after eq. ([4.31|) , the residues Zs{p) are positive functions. 

Let us construct now covariant gauge polarization vectors appropriate for the lon- 
gitudinal (e^(p ; 0)), and, respectively, for the transverse (e^(p ; A), A = 1, 2) soft normal 
modes. They satisfy 

*D,J(.;,(p),p)e^(p;A) = 0, (A.7) 

where uj\{p) is equal to uji{p) for A = and, respectively, to uJt{p) for A = 1, 2. By 

writting e'^(p; A) = (0,e(p; A)), it follows that P^e^(p; A) = 0, for A = 0,1,2 and p° = 

uj\{jp)] furthermore, p ■ e(p; A) = 0, for A = 1, 2, and {5^^ — p''p^)e^{p; 0) = 0. With the 
normalization e(p; A) ■ e*(p; A') = — 6xx', we obtain 

e^(p;0) = (3{p) (1, pooi/p), /3(p)/3*(p) = p'/ioof-p'), (A.8) 

and 

e'^(p; A) = (0, e(p; A)), p ■ e (p; A) = 0, e (p; A) ■ e* (p; A') = 5aa', (A.9) 

for A = 1, 2. Note also that 

Y: e'^lp; A) e^*{p; A) = Vf , e^{p; 0) e^*{p- 0) = Vf ■ (A.IO) 

A=l,2 
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The soft gluon effective propagator ( |A.3| ) is intimately related to the dielectric per- 
mitivity tensor, commonly used in classical plasma physics to describe the polarization 
properties of the plasma. In the weak field (or the abelian) limit, where the plasma re- 
sponds linearly, (i.e., j"^'^ = Yl^^A'^), we can use the standard definition of e*-']!^. The 
polarization vector P{u!, p) satisfies, by definition, —iujP{u!, p) = j*"°'(co', p), and also 
zp ■ P{uj, p) = — p*"°'(u;, p), because of the continuity equation for j^"°', eq. (|2.10|) . The 
linear relation between the displacement vector D(ci;, p) = E(c<j, p) + P(u;, p) and the 
electric field defines the dielectric tensor: 



D\co,p)^^^ico,p)E^icu,p). 



(A.ll) 



It follows immediately that 



e'^{u;,p) = 6'^ --W^iuj,p). 



u^ 



(A.12) 



For isotropic dielectric media, e*-^ has only two independent components, and it is useful 
to choose them as the longitudinal and the transverse dielectric functions, defined by 

(A. 13) 



e'^(tu, p) = ei{uj,p)py +et{uj,p) [5'^ - py ^ 
From this relation, together with eqs. ( |A.12| ) and ( p.l2| ), it follows that 



e/(^,p) 



1 - 



uj"^ — p2 



nK^,p) 






Qiu/p) 



(A. 14) 



1 / N 3 cj^ 

et{uj,p) = 1 -Ut{uj,p) = I - - — 

a;^ 2 p-^ 



where 



Q(a;) = ^ln^±l = - (In 
^^ ' 2 x-l 2 



x+ 1 



X — 1 



1 - ^^^ Q{u;/p) 



Z7r^(l-|x|)') 



(A. 15) 



(A.16) 



and the prescription x -^ x + irj, t] —>■ 0^, has been understood in writting the imaginary 
part of Q{x). Remark the imaginary parts in ei^t{^,p), which occur for \uj\ < p, as 
expected. From eqs. (|A.14|) , ( [A.15|) and ( |4.30| ), we derive 



^(..Q(-,P)) 



^i \p) duo 



H 



^i (p) ujf (p) — p^ 



- 1, 



(A. 17) 



and, similarly, 

d 

dui 



uj{et{uj,p) -p'^/uj^) 



^t{p) 



^t{p) -7- {et{^,p) -pV^^) 



duo 

H 



^tip) 



+ 



p^ 



^tip)-p^ ^t{p) 



1. 



(A.18) 
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The dielectric functions are related to the inverse effective propagators, as shown by 
eqs. O and (|X1^ )- (|A:T5D , 



* A -1/ 



^l {^,P) = {^ -P)^1{^,P), *^t {^,P)=^ ^t{^,P)-P ■ 



(A.19) 



From the defining relation of the residues, eq. ( |A.6| ), it follows that 



^lip) 



dei 

duo 



^ 2uf{p) 
^lip) ojfip) — P' 



^r\p), 



u;t{p) T- (6t(^,p) -pV^') L . . = 2^r'(p)- 



^t(p) 



(A.20) 



Note finally the following angular integrals, which are needed in Sec. 4. 3 in order to 
compute the energy and the momentum of the plasma waves: 



dQ v' v^ 



and 



An (to — V ■ p)2 



dQ v^ v^ 

An (cj — V ■ p) 



aiuj/p) : : b(^/p) , rin -i -1\ 



p^ 



p^ 



^ 3 PV + 3 {5'^ -p'p>). 



po 



po 



(A.21) 



(A.22) 



These results are valid for uo > p. In these formulae 



aix) = 1 + 



X 



x'^ — 1 



2Q(x), h{x)=Q{x) - 1, 



and 



X 



c{x)^-\q{x) - (l-^ 

X \ \ a;^ — 1/x"^ — 1 



1 / x^ 
d{x) = — 



2x \ a;^ — 1 



(A.23) 



Q{x) , (A.24) 



with Q{x) given by ( [A.16| ). For a; > 1, the functions a(x), 6(x), c(x) and d{x) are all 
positive and vanish only in the limit x -^ oo. 
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